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Abstract 

We have studied the production of B hadrons in 1.8-TeV pp cohisions. We 
present measurements of the fragmentation fractions, /„, f^, fs and fbaryam 
of produced b quarks that yield B^, B^ and hadrons. Reconstruction 
of several electron-charm final states yields fs/{fu + fd) = 0.213 ± 0.068 and 
fbaryon/{fu+fd) = 0.118±0.042, assuming /„ = fd. If all B hadrons produced 
in pp collisions cascade to one of these four hadrons, we determine fu = fd = 
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0.375 ± 0.023, /, = 0.160 ± 0.044 and fbaryon = 0.090 ± 0.029. If we do not 

assume = fa, we find fa/fu = 0.84 ± 0.16. 

PACS Numbers: 13.60.Le, 13.60.Rj, 14.65.Fy, 13.87.Fh 
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Typeset using REVT^iX 



Bottom (6) quarks are not observed as independent entities but are confined with a 
partner antiquark or diquark inside hadrons. Once a b quark is produced, the process by 
which it combines with quarks and gluons to form a hadron is called fragmentation and 
is governed by the strong force, described by the theory of Quantum Chromodynamics 
(QCD) ||T]]. In this fragmentation process, the color force field creates additional quark- 
antiquark partners that then combine with the bottom quark to create a B hadron. 

The process by which a bottom quark fragments into a hadron cannot be reliably cal- 
culated using perturbative QCD methods. Therefore, the fragmentation properties of the 
b quark must be determined empirically. In this Letter, we investigate one such property, 
namely the flavor dependence of the fragmentation process for bottom quarks produced 
in 1.8-TeV pp collisions. Our results provide the most accurate measurements of this fla- 
vor dependence and for the first time bring together in one study all previously studied B 
hadrons. 

We define fu, fd, fs and fbaryon to be the probabilities that the fragmentation of a 6 
quark will result in a weakly decaying _B+, 5°, 5° meson and A° baryon, respectively. We 
explicitly include in these "fragmentation fractions" contributions from production of heav- 
ier B hadrons that decay into final states containing a B^, B^, B^ meson or A° baryon. 
The ALEPH experiment used reconstructed B^ D^l^vX decays produced in e^e" colli- 
sions at the Z° resonance to determine the value fs = 0.120to;o34 The LEP Working 
Group on B Oscillations has compiled B^^ mixing results from the four LEP experiments 
and the Collider Detector at Fermilab (CDF) experiment for the mixing parameters x and 
Arriii 0. The average values of these parameters constrain the value of fs, yielding the result 
/, = 0.10lt[j:|]?[i |[. The CDF experiment has measured + fd) = 0.210 ± 0.036tE];|]i 
using double semileptonic decays b cfiX with c — > sfiX |^. The ALEPH and DELPHI 
experiments measured fbaryon by reconstructing A° A'l^uX decays PJ§]. Their combined 
result is fbaryon = 0.10lto;o3i 0. The CLEO experiment determined the quantity analogous 
to fd/fu, fVf^ = {T{4S) B^B^)/{T{4S) B+B') = 0.88 ± 0.16 and 0.90 ± 0.14, by 
reconstructing B D*liJ decays 0] and B J/il)K'^*'> decays 0], respectively. Both of 
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these measurements are consistent with the isospin symmetry expectation that fd = fu- 

Our measurement is performed by reconstructing B hadron semileptonic decays to elec- 
trons and charm hadrons from a 107 pb^^ sample of 1.8-TeV pp collisions recorded by 
CDF during 1992-95. The ratios of the b quark fragmentation fractions, namely fd/fu, 
fs/ifu + fd) and fbaryon/ifu + fd), are determined from the B hadron production ra- 
tios. We reconstruct the B hadrons in the following decay modes and their charge con- 
jugates: B+ Tfe+UeX where D° ^ K+n-; 5° D*'e+UeX where D*- ifn- and 
D° ^ K+n-; B° D-e+UeX where K+n-n-; ^ Dje+UeX where D; ^vr" 

and (f) — > K~^K~] and A° — * A^e^u^X where A~ — > pK^n^. The average transverse 
momentum of the B hadrons we reconstruct is 20 GeV/c. 

The CDF detector has been described in detail elsewhere The CDF coordinate sys- 
tem defines the z axis along the proton beam direction and the polar angle 6 with respect 
to the z axis. The azimuthal angle is measured in the plane transverse to the beam. 
The relevant detector components for this measurement are the charged-particle tracking 
system and the calorimeters. The tracking detectors lie inside a 1.4-T solenoidal magnetic 
field. The silicon vertex detector (SVX), positioned immediately outside the beampipe, pro- 
vides precise charged particle reconstruction and allows identification of displaced vertices 
from secondary decays. The central tracking chamber (CTC), which encompasses the SVX, 
measures the momenta of charged particles over a pseudorapidity range \7]\ < 1.1, where 
rj = — lntan(^^/2). The central electromagnetic (CEM) and hadronic (CHA) calorimeters, 
arranged in a projective tower geometry, surround the tracking volume and are used to mea- 
sure clusters of energy deposited by electrons, photons and hadrons. The central electro- 
magnetic strip chamber (CES), embedded in the CEM at the position of shower maximum, 
measures the electromagnetic shower profiles in the and z directions. 

A three level trigger system is used to identify electron candidates, with the first level 
requiring a CEM energy deposition greater than 8 GeV. The electron candidates sat- 
isfy a Level 2 trigger that requires a spatial match between a track in the CTC with 
Pt > 7.5 GeV/ c and an energy cluster in the CEM with Et > 8.0 GeV, where Pt = PsinO 
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and Et = EsinO. The fraction of hadronic energy in the cluster is required to be smalL 
We require a spatial match of the CTC track to a cluster of energy in the CES and apply 
a threshold requirement to the energy deposition in the CES. The Level 3 trigger requires 
that the lateral shower profiles in the CES and CEM be consistent with those expected for 
an electron, and re-applies the previous trigger criteria with greater precision. Approxi- 
mately six million electron candidates survive this trigger selection. We reduce the sample 



to three million candidates by applying more stringent criteria [|10l- We require that the 
fraction of hadronic energy in the cluster be less than 4%. We reject electron candidates 
that are likely to be from photon conversions and from and boson decays. Finally, to 
ensure a uniform electron identification efficiency in the different B hadron decay topologies, 
we reject candidates with more than one track pointing at the CEM cluster and demand 
that the ratio of cluster energy to track momentum be in the range 0.75 < Et/ Pt < 1.40. 

The semileptonic B hadron decays are identified by reconstructing the charm hadron in 
the vicinity of the electron. The meson is reconstructed by identifying the products of 



the fT+TT" decay in a cone R = ^ (Ar^)^ + (A0)^ = 1.0 around the electron track. The 

charge correlation between the electron and the charm hadron daughters from semileptonic B 
hadron decays is exploited to reduce the contamination from random combinations of leptons 
and charmed hadrons. Particles arising from the weak decay of a i? hadron are normally 
displaced from the primary vertex. Therefore, we require the charm-hadron daughter tracks 
to have an impact parameter [do) inconsistent with zero {\do\/a{do) > 1.5, where a{do) 
is the uncertainty on do). The combinatorial background is further reduced by requiring 
that Pt{K) > 1.2 GeV/ c and Pt{t^) > 0.5 GeV/ c, which are the same criteria used in the 
reconstruction of the other channels, except where noted. The daughter tracks are required 
to be consistent with coming from a secondary vertex that is displaced in the transverse 
plane from the pp interaction point {L^y/ a{Lxy) > 1). Finally, the invariant mass of the 
electron- charm system is required to be less than 5.0 GeV/c^. 

The invariant mass of the Kn candidates in the electron sample is shown in Fig. |I|(a). 
To this distribution we fit the sum of a Gaussian signal and an exponential background and 
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count 1848 ±58 D signal events. We observe no significant signal in the combinations with 
the wrong electron-hadron charge correlation. 

The D*~ meson is reconstructed in the D*^ D°7r~ channel. We consider candidates 
with 1.80 < M{K7t) < 1.95 GeV/ c^, where M is the mass, and consider all charged particles 
with Pt > 0.4 GeV/ c for the additional pion. The mass difference distribution, AM = 
M{K'K'k) — M{Kit), is shown in Fig. 0(b). To this distribution we fit a double-Gaussian 
signal and a background modeled by a threshold function. We reconstruct 249 ± 19 D*~ 
signal events. 

The D~ meson is reconstructed in the K^h^tt" channel. In this channel, the 

three daughter tracks are required to form a vertex. The invariant mass distribution of the 
Ktiti candidates is shown in Fig. 0(c). To this distribution we fit the sum of a Gaussian 
signal and a linear background and count 736 ± 62 signal events. 

The meson candidates are identified by looking for the products of the D~ 07r~ 
decay, where — K^K~ . Both kaons and the pion are required to come from a common 
vertex. This decay chain provides two additional criteria effective in rejecting combinatorial 
backgrounds. First, we require that the mass of the K" system be within 0.010 GeV/ (? 
of the world average (f) mass of 1.019 GeV/ (?. Second, we impose the criterion | cos^/^l > 0.4, 
where is the helicity angle between the and K meson candidates in the rest frame. 
The invariant mass distribution of the KKtt candidates is shown in Fig. [l|(d). We reconstruct 
59 ± 10 signal events. 

The baryon candidates are identified by looking for the products of the A~ pK'^n' 
decay. We require Pt{p) > 2.0 GeV/c. Since the relative combinatorial background under 
the A~ signal is large, we also require that the specific ionization (dE/dx) deposited by the 
proton candidate in the CTC be consistent with that expected for a proton. The invariant 
mass distribution of the pK it candidates is shown in Fig. |l](e). We reconstruct 79 ± 17 A~ 
signal events. 

The D~e^ and A^e+ final states represent relatively pure samples of the i?^ and A° 
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hadrons, respectively. However, the remaining electron-charm final states that we recon- 
struct come from several i?-meson species through feed-down from vector and orbitally- 
excited D meson decays. For example, the decay 5^ Dl*'e^Ve can be followed by the 
decay D**" D~K^. This channel contributes to the D^e^ sample but refiects produc- 
tion rather than production. We use the branching fractions for each decay to determine 
the feed-down contributions. These branching fractions are derived from the measured val- 
ues according to the spectator model and isospin symmetries [0 . 

The spectator model predicts that the inclusive semileptonic decay widths for the various 
B hadrons are equal, yielding, for example, the relation 

{nB)) 

where B represents the branching fraction and r is the lifetime. A similar relationship holds 
for the exclusive semileptonic branching fractions for the three B mesons. We use isospin 
symmetry to calculate the branching fractions for the D* and D** decays that feed down 
into the final state D mesons that we reconstruct. 

The acceptance and reconstruction efficiency for each final state vary according to 
whether the D meson is produced directly in the B meson decay or is the daughter of 
an excited D meson state. Several efficiencies, such as the electron identification efficiency, 
the conversion removal efficiency and the two-track finding efficiency, cancel in the ratio 
of fragmentation fractions. Of the remaining efficiencies, the single-track finding efficiency 
and the electron trigger efficiency are measured using CDF data. We use a Monte Carlo 
calculation to determine all other acceptances and efficiencies. This uses a next-to-leading 
order perturbative calculation of the differential cross section for 6-quark production in pp 
collisions [|11|] followed by fragmentation governed by the Peterson formulation We use 
the ISGW |13| model to determine the semileptonic B hadron decay kinematics ||14|| . 



The systematic uncertainties on the reconstruction efficiencies include those associated 
with the tracking and trigger efficiencies and Monte Carlo statistics. We also assign an 
uncertainty to account for the poor knowledge of the A|| production polarization in pp 
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collisions. Finally, we consider the possibility that the Peterson fragmentation parameter eb 
may differ for each B hadron species. We assign the fractional uncertainties of 6.4% and 
6.1% to account for the possible variation of eb for and A° hadrons, respectively. We 
determine these values by calculating the larger variation in the reconstruction efficiency for 
values of one standard deviation away from a central value of eb = 0.006 ± 0.002 p!5| . 
These contributions represent the largest uncertainties associated with the reconstruction 
efficiencies. The total fractional systematic uncertainties are ±2.1%, ±3.4%, ±5.4%, ±7.3%, 
and ±8.7% for B+ Tfe+u^X, 5° ^ D'e+UeX, 5° ^ L)*-e+z/eX, 5° ^ D^e+UeX and 
A° A^e^UgX decays, respectively. 

In order to determine the fragmentation fractions taking into account the cross contam- 
ination and feed-down, we fit the five observed event yields and their uncertainties to the 
three ratios of fragmentation fractions. We formulate the problem by defining a function 
comparing the predicted with observed event yields. We allow the semileptonic branching 
fractions for the B mesons to vary in the fit, constrained to their measured or calculated un- 
certainties. We note that the measured branching fractions include the implicit assumption 
that /7/+ = 1. 

We make this measurement assuming isospin symmetry by fixing fd/ fu = 1 in the fit. The 
fit results in the values /,/(/„ ± fd) = 0.213 ± 0.068 and haryon/Uu + fd) = 0.118 ± 0.042, 
where uncertainties on the event yields, the reconstruction efficiencies and the branching 
fractions are included. We can determine the absolute fragmentation fraction values from 
our fits by assuming that the B^ , B^, B'^ and hadrons saturate the 6-quark production 
rate, i.e., fu + fd + fs + fbaryon = 1- This relationship yields fu = fd = 0.375 ± 0.023, 
fs = 0.160 ± 0.044 and fbaryon = 0.090 ± 0.029. By incorporating another term in the 
function, we have combined these results together with the complementary measurement 
by CDF of fs using double semimuonic decays to determine the more precise values 
of f^ = fa = 0.375 ± 0.015, /, = 0.160 ± 0.025 and fbaryon = 0.090 ± 0.028. Results 
for the fragmentation fractions obtained using all available measured exclusive semileptonic 
branching fractions instead of employing the spectator model predictions are consistent with 
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the results presented here. 

By relaxing the isospin symmetry requirement, we find that fd/fu = 0.84 ± 0.16, con- 
sistent with isospin symmetry and the measurements of /°/ by the CLEO collaboration. 
The individual values for fg and fbaryon remain unchanged. 

In conclusion, we have measured the four b quark fragmentation fractions for weakly 
decaying B hadrons produced in pp collisions. We measure fbaryon = 0.090 ± 0.029, in good 
agreement with measurements made on B hadrons produced in high energy e'^e" collisions 
at LEP. The pp result, fg — 0.160±0.025, is two standard deviations higher than the current 
world average value, which is dominated by LEP measurements and by inference from B^^ 
mixing measurements. 
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FIG. 1. Invariant mass distributions of charm hadron candidates in 107 pb^ of inclusive elec- 
tron data, a) Ktt invariant mass distribution for candidates, b) Mass difference distribution, 
AM = M^Kirir) — M{Kit), for D*^ candidates, c) Ktttt invariant mass distribution for D~ 
candidates, d) KKtt invariant mass distribution for candidates, e) pKn invariant mass dis- 
tribution for candidates. The shaded histograms represent the combinations with the wrong 
electron-hadron charge correlation. The shaded area in a) has been scaled by 0.5 for display 
purposes. 
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